FIRBACK FAR INFRARED SURVEY WITH ISO: 
DATA REDUCTION, ANALYSIS AND FIRST RESULTS 
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ABSTRACT 



FIRBACK is one of the deepest cosmological surveys 
performed in the far infrared, using ISOPHOT. We 
describe this survey, its data reduction and analysis. 
We present the maps of fields at 175 fim. We point 
out some first results: source identifications with ra- 
dio and mid infrared, and source counts at 175 fim. 
These two results suggest that half of the FIRBACK 
sources are probably at redshifts greater than 1. We 
also present briefly the large follow-up program. 



of today in the field of galaxy formation and evolu- 
tion (see Puget & Lagache, this volume). 
In this context, FIRBACK, which is a cosmological 
survey dedicated to a study of the Cosmic Far In- 
frared Background and galaxies contributing to it, is 
one important step. 



2. DATA SET 



Key words: ISO; ISOPHOT; far infrared survey; cos- 
mology; galaxy evolution 



1. INTRODUCTION 



FIRBACK is a survey of 3 main fields (described in 
Table [j]) covering 4 square degrees using ISOPHOT. 
All fields have been observed in raster mode (AOT 
P22) with the C200 detector and C_160 filter centered 
at A = 175 \xm. We have also additional ISO obser- 
vations with ISOPHOT and ISOCAM (see below). 



The discovery of the Cosmic Far Infrared Background 
Radiation fCFIBR 
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Lagache et al 
COBE data and constraints on its spectrum 
strates that about two thirds of the light emitted by 
galaxies integrated over all redshifts has been pro- 
cessed by dust and released at FIR and submm wave- 
lengths. This background is in line with models using 
strong evolution of galaxies. 

Studying the population of galaxies radiating mostly 
at long wavelengths as a function of redshift is thus 
one of the important but difficult observational tasks 



Table 1. Fields of the FIRBACK survey. 
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Figure 1, Map of FIRBACK Marano fields 
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Figure 2. Map of FIRBACK ELAIS Nl field 



FIRBACK Marano 

This area (map Figure 1) is composed of four indi- 
vidual fields, called Marano 1, 2, 3 and 4, which have 
been observed four times. The redundancy goes from 
32 (in overlaping regions) to 1 (on the edges) with 
an average of 16, and the oversampling is optimal 
in both Y and Z directions in Marano 2, 3 and 4 



(half pixel offset between rasters). Marano 1 covers 
31.5' x 31.5', and Marano 2, 3 and 4 about 77' x 26'. 
Integration time per sky pixel is 256 s on average. Ad- 
ditional observations: in Marano 1 we have one P25 
absolute measurement, five small rasters at 90 jim 
(not yet reduced); in the complete Marano field, we 
have a 15 fim ISO CAM survey. 



FIRBACK ELAIS Nl 



3.2. Map Projection 



This area (map Figure 2) is composed of eleven indi- 
vidual fields, observed 2 times with an offset of less 
than a pixel, and the redundancy goes from 16 (in 
overlaping regions) to 1 (on the edges) with an aver- 
age of 8. ELAIS Nl covers about 2 square degrees. 
Integration time per sky pixel is 128 s on average. 
Note that the northern FIRBACK fields have been 
choosen to cover some corresponding ELAIS fields at 
other wavelengths (see Rowan- Robinson et al., this 
volume) . 



FIRBACK ELAIS N2 

This area (map Figure 3) is composed of nine indi- 
vidual fields, observed 2 times with an offset of one 
pixel, so the redundancy is 8 on average, without 
oversampling. ELAIS N2 covers about 57' x 57'. 



Our basic data projection is performed in two steps: 
(1) geometrical grid computation and (2) projection 
of the signal. 

- The first step is the most time consuming. We read 
all AAP data of one complete field, and create a large 
grid in (a, S) with tunable pixel sizes (usually 10 arc- 
sec); this step has to be performed just one time. 

- The second step just acts to restore the grid, to 
coadd individual raster data onto the grid by dis- 
tributing the pixel signal in the area of one pixel, 
and by updating a weighting map. 

This method is well-suited to large amounts of data, 
in terms of sky surface or number of co-added rasters 
with different roll angles. Figures |], || and |3| show 
such maps. 



4. RESULTS 



3. DATA REDUCTION 



4.1. Source Detection 



3.1. From Raw Data to AAP 



We used PIAR version 7 .2.2 (Phot Interactive Anal- 
ysis, pabriel et al. 1997 ) to reduce data from ERD 
(Edited Raw Data, the lowest level) to AAP (Auto 
Analysis Product Data, level where flux density and 
positions are available): 

- ERD to SRD (Signal per Ramp Data): linearisa- 
tion, deglitching, fit ramps to first order 

- SRD to SCP (Signal per Chopper Plateau Data): 
deglitching, no drift correction 

- SCP to AAP: no reset interval correction (reset 
times are identical), subtraction of dark current, 
power calibration using interpolation between 2 FCS. 
We then used our own IDL routines to correct for 
transients induced by strong glitches, and to perform 
a flat correction using observational redundancies. 



X PIA is a joint development by the ESA Astrophysics Di- 
vision and the ISOPHOT Consortium 




Source detection is directly performed on final co- 
added maps, using eye recognition to discriminate 
point sources from cirrus structure. We then com- 
puted the flux by aperture photometry with radii 
choosen to minimize noise in source detection R m in — 
60" and R max = 92" flPuget et al. 1999Q. 



4.2. Catalogs and Identifications 

Sources in our catalog have a S/N ratio greater than 
3, and photometric accuracy is about 30 %. Fluxes 
are in the range 75 mJy to 1.6 Jy, with an average of 
252 mJy and a median of 175 mJy. 
In order to constrain both Spectral Energy Dis- 
tribution (SED) and position of these sources, we 
looked for identifications: Table E gives a summary 
of detected sources and identifications with ISO or 
ground-based radio telescopes. 

In the Marano Field, we performed a smalle r 15 /jm 
survey using ISOCAM ( besert et al. 199SJ ), which 
covers about 50 % of the 175 \im survey surface. 
Radio identifications are performed usin g Australia 



Telescope Compact Array available data (Gruppioni 



at al. 1997 ) or preliminary data of our program being 
reduced, that is why 8 identifications is a lower limit. 
The situation is different in northern fields ELAIS 
Nl & N2 because the ELAIS consortium performed 
observations using ISO at A = 15 /xm, A = 90 /im 
(Rowan-Robi nson et al.. this volu me), using the VLA 
at A = 20 cm ( Ciliegi et al. 1998 ), and at other wave- 
lengths. Most of the 31 sources with 90 /im coun- 

-F l75n 



terpart have a color ratio 



> 2, which 



^90 ptm PHOT 

has to be compared to local sources (z < 0.1) where 
Fl75 " m 1 dStickel et al. 1998|), suggesting that 



^lOO jj.m IRAS 

a significant proportion of FIRBACK sources are 
large-redshift sources (see discussi on below) . Most o f 
NEDR identifications are FIRST ( [Becker et al. 1995| ) 
or IftAS sources. 



Figure 3. Map of FIRBACK ELAIS N2 field 



2 The NASA/IPAC Extragalactic Database (NED) is oper- 
ated by the Jet Propulsion Laboratory, California Institute of 
Technology, under contract with the National Aeronautics and 
Space Administration 
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Table 2. Number of FIRBACK sources and identifi- 
cations 



Field Name 


175 fim 


15 fim 


90 /im 


20 cm 


Fu(mJy) > 


75 


1 


60 


0.2 


Marano 
ELAIS Nl 
ELAIS N2 


78 
113 
85 


29 

19 

no data 


no data 
18 
13 


> 8 
20 
22 



We 



4.3. Source Counts 



present raw source counts in Figure 0, using 
sources coming from FIRBACK ELAIS Nl. !~N2 and 
Marano 2 3 4 fields (208 sources brighter than 100 
mjy in less than 3 square degrees). Raw counts 
means here that counts are not corrected for incom- 
pleteness. For log(S > 300 mjy), counts are not ex- 
pected to change a lot due to incompleteness correc- 
tion, so the steep slope of 2.20 is a minimum. We 
also plotted Marano 1 raw counts (Puget et al. 1999) 
which are higher by a factor 2 but compatible witmn 
uncertainties, probably because based on 22 sources 
in a smaller area (0.25 square degree), Kawara et al. 
(1998) data in the Lockman Hole (difficult to see in 
the plot because their count at 150 mjy gives the 
same number as ours), models with and without evo- 
lution from Franceschini et al. (1998), and model E 
from Guiderdoni et al. (1998). Compared to extrap- 
olated IRAS counts, FIRBACK counts are about 20 
times higher at 200 mjy. Our counts are not com- 
patible with non-evolution models and need strong 
evolution. Model with evolution from Franceschini 
et al. (1998) predicts smoother slopes than observed, 
whereas model E from Guiderdoni et al. (1998) fits 
roughly the slope even if lower than observed. 



5. FOLLOW-UP | 

An intensive multiwavelength follow-up program is 
being performed. 

In the Marano field, we are following-up at A = 20 
cm using the Australia Telescope Compact Array in 
compact configurations (data being reduced) and ex- 
tended configuration (program to submit) (Dole). In 
the optical range, one source has been observed with 
the NTT (Dennefeld). 

In the ELAIS Nl & N2 fields at millimetre wave- 
lengths, we are following-up using the IRAM 30 m 
antenna and Plateau de Burc Interferometer (Omont 
& Guiderdoni). Submillimetre observations in Nl are 
scheduled at JCMT (Scott) and in April at CSO (La- 
gache) , and optical follow-up of sources is performed 
at Palomar (Reach). 



DISCUSSION ON COSMOLOGICAL 
IMPLICATIONS 




3 extended FIRBACK collaboration for follow-up also in- 
cludes: D. Benford, P. Cox, M. Dennefeld, G. Helou, R. 
Mac Mahon, A. Omont, F. Pajot, T. Phillips, D. Scott, F. 
Viallefond 
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Figure 4- Source counts with FIRBACK data in 
Marano 1, 2, 3, 4, ELAIS Nl and ELAIS N2. These 
counts are not corrected for incompleteness. Data: 
diamonds: Marano 1 (Puget et al. 1999); triangles: 
ELAIS Nl, N2 and Marano 2 3 4; Lockman Hole data 
(Kawara et al. 1998): purple star at 150 mJy, su- 
perimposed on our counts. Models: Guiderdoni et al. 
1998: blue line: model E; Franceschini et al. 1998: 
red dash-dash: with evolution; green dash-dot: with- 
out evolution 



6.1. Population 

The steep slope of the number counts (> 2.2) cannot 
be accounted for by the effect of the K-correction (ra- 
tio, at a given wavelength, of the emitted flux and the 
observed flux with a spectral redshift) if no cosmo- 
logical evolution is present as can be seen in Figure 0. 
The Guiderdoni et al. (1998) model which gives 
rather a good fit to our counts presents a strong evo- 
lution for the star burst component (oc (1 + z) 5 ) and 
an ULIRG (Ultra Luminous Infrared Galaxies) com- 
ponent. In this model, the far infrared background 
comoving luminosity per unit volume increases by a 
factor of 20 from z — to z = 2.5. 
This model also predicts a redshift distribution for 
sources brighter than 100 mjy at 175 fim. The red- 
shift distribution peaks at z = 0.9 with a median of 
z = 1.0, about half of the sources lie at redshift be- 
tween 1 and 2, and about 10% above z = 2. Further- 
more, this model is compatible with the far infrared 
backgound detected in the COBE data. Sources de- 
tected in our survey account for 3% of this back- 
ground at 175 fim, and the model predicts that 90% 
of the background is made of sources brighter than 
2 mjy; this predicted value is probably too low, be- 
cause our slope is higher than the model, and our 
counts are not yet corrected for incompleteness. 

6.2. Identifications 
If these galaxies are star burst galaxies with a SED 
peaking at 75 fim, a typical source observed close to 
our detection limit (£175 Mm = 100 mjy) has a lu- 
minosity 1.2 x 10 12 L Q (in a h = 0.5 and fio = 1 
cosmology) . 
The same source will have P15 ^ m ~ 1 mjy and thus 




Figure 5. K-correction at 175fj,m (black line) and 
90 fim (blue dash-dash line) as a function of redshift 
for a starburst galaxie whose SED peaks at 75 fim 



is easily detectable with ISOCAM. One of the differ- 
ences between mid- and far-infrared observations is 
that the K-correction allows observation of galaxies 
with redshifts in general below 1.4 (because of the 
PAH feature's cutoff) in the mid infrared; no such 
strong limit exists in the far infrared. We have only 
25% detections at 15 fim, implying that more than 
half of the sources are at z > 1.4 (redshift at which 
15 fim flux drops sharply). 

The expected 90 fim flux for this source is -Fgo Mm — 
40 mJy, and thus not easily detectable in the ELAIS 
survey. Futhermore, the K-correction at 90 and 
175 fim as function of redshift (Figure ph shows that 
between redshifts z = 1 and z = 2 the ratio 
-j7 ' grows from 6 to 14: at z — 2 a simi- 

lar galaxy has 14 times more flux at 175 fim than at 
90 fim. 

The non-thermal emission is known to be strongly 
correlat ed with the FIR e mission for starburst galax- 
ies (e.g |Helou et al. 1985 ), even if Sanders & Mirabel 
(1996) show an important scatter in the correlation. 
With ELAIS VLA identifications (42 sources), we 
compute q = log{F vFIR ) - log(F„i aghz) (we take 
Firs fim for F v fir and not a combination of Fqq ^m 
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and F WOtim ) the FIR to 1.4 GHz log ratio. For 70 % 
of our sample q < 2.3, and 2.3 < q < 3 for 30 %; 
the dispersion is significant and must be taken into 
account for deep radio observations. 



7. CONCLUSION 



We presented the FIRBACK 175 fim survey, data re- 
duction and analysis. Identifications of sources at 
various wavelengths, color ratios and source counts 
suggest that half of the sources are at redshifts 
greater than 1. Following-up these sources in order 
to have better positional accuracies, and SED con- 
straints is now the key for understanding this popu- 
lation of dust enshrouded galaxies. 
Additional information, images and papers are avail- 
able online at http://wwwfirback.ias.fr 
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